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I.  INTRODUCTION 

For  several  years  theoretical  models  (References  1 and  2),  which  predict 
the  bound-bound  infrared  (1-10/tml  emissivity  of  an  atmospheric -type  plasma, 
have  been  used  with  little  experimental  comparison.  These  calculations,  done  at 
Naval  Research  Laboratory  (NRL ) and  other  places,  use  a collisionaJ-radiative 
model  (References  3 through  7)  to  calculate  the  number  densities  of  the  excited 
states  of  oxygen  and  nitrogen  atoms  submerged  in  a plasma.  From  the  number 
densities  of  the  excited  states  of  an  atom,  one  can  easily  calculate  the  bound- 
bound  optically  thin  emission  in  any  spectral  region  from  (Reference  8) 


I = — A N L . 
4 7T  Ul  U U 


(1) 


where  I is  the  intensity  of  the  radiation,  u.’  is  the  angular  frequency  of  tiie  radia- 
tion emitted.  AuJ  is  the  atomic  transition  probability,  Nu  is  the  number  density 
of  the  excited  state,  ;uid  L is  the  length  of  the  plasma. 

The  purpose  of  this  experiment  is  to  provide  a check  of  the  various  assump- 
tions which  go  into  the  theoretical  models  mentioned  above.  To  do  this,  the 
number  densities  of  selected  excited  states  of  oxygen  and  nitrogen  atoms  sub- 
merged in  a non-equilibrium  plasma  have  been  measured.  This  was  done  by 
measuring  the  absolute  intensity  of  bound-bound  visible  radiation  emitted  by 
atoms  undergoing  transitions  from  the  state  of  interest  to  some  other  state 
(References  9 through  11 ) . This  combined  with  good  plasma  diagnostics  to 
measure  the  plasma  length  is  then  used  along  with  equation  (1 ) to  give  the 
excited  state  number  densities.  Of  course  the  spatial  development  of  the  plasma 
parameters,  electron  density  and  temperature,  must  also  be  known  so  that  mi 
unambiguous  comparison  of  the  experimental  number  densities  can  be  made  with 
those  generated  by  the  computer  programs. 

The  atmospheric  simulation  computer  programs  were  originally  developed 
to  investigate  the  collisional-radiative  recombination  of  a low  electron  density 

0103 

(Np  ~ 10  - 10  cm"  ) plasma  with  a temperature  of  a few  tenths  of  an  electron 
volt.  In  practice,  it  is  difficult  to  detect  the  line  radiation  emitted  bv  a plasma 
at  these  low  densities  and  temperatures  when  the  size  of  the  plasma  is  restricted 
to  laboratory  dimensions.  For  this  reason,  our  measurements  were  performed 
at  electron  densities  between  1012  and  1014  cm"3  with  temperatures  between  0.3 
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and  1.5  eV.  It  should  be  noted  that  the  assumptions  made  when  developing  the 
model  codes  should  still  be  valid  in  this  region  thereby  making  a comparison 
possible. 

The  next  section  will  describe  the  () -pinch  used  as  a plasma  source  in  this 
experiment.  After  that  the  oxygen  diagnostics  and  results  will  be  discussed  in 
total.  Then  the  nitrogen  diagnostics  and  results  will  be  presented.  The  final 
section  will  be  a summary  of  results  and  some  comparisons  between  the  two 
experiments . 

II.  APPARATUS 

The  plasma  source  used  for  this  experiment  was  a (> -pinch.  A hybrid  sche- 
matic is  shown  in  Figure  1 . The  tube  is  made  of  quartz  and  is  137  cm  long  with 
an  inside  diameter  of  9 cm.  The  single-turn  d-pinch  coil  is  80  cm  long  and  has 
a slot  in  the  center,  for  side-on  viewing,  which  is  not  shown  in  the  diagram.  The 
ambient  gas  in  the  discharge  tube  is  preionized  by  discharging  a free  ringing 
0.126  fi F.  20  kV  capacitor  through  small  side  coils.  The  plasma  is  then  further 
heated  by  discharging  a low  inductance  free  ringing  1 .5  j.  F,  20  kV  capacitor 
through  the  single-turn  main  coil.  Finally  the  high  inductance  250  pF.  20  kV 
main  bank  is  discharged  through  the  single-turn  main  coil.  This  bank  is  clamped 
and  has  a quarter  period  of  about  18  psec  with  an  e-folding  decay  time  of  about 

200  psec.  In  this  manner,  a slowlv  decaving  plasma  is  produced  with  initial 

15  -3 

conditions  of  electron  density  ^10  cm  and  electron  temperature  "3  eV.  The 
large  initial  field  - 9 kG  produced  by  the  main  bank  discharge  helps  isolate  the 
walls  of  the  quartz  tube  from  the  hot  early  plasma.  Thi.:  helps  to  improve  the 
purity  of  the  plasma.  The  longevity  of  the  magnetic  field  produced  by  the  main 
bank  discharge  improves  the  radial  homogeneity  of  the  plasma  by  slowing  down  the 
radial  loss  rate. 

III.  OXYGEN  MEASUREMENTS 

A.  Diagnostics 

The  ambient  fill  tins  used  throughout  this  part  of  the  experiment  was  pure 
O2  at  a pressure  of  10  mTorr.  The  first  diagnostic  done  on  the  plasma  was  to 
take  the  time-integrated  spectrum  shown  in  Figure  2.  The  lower  spectrum  is  that 
of  the  plasma  emission  observed  down  the  axis  of  the  0 -pinch  tube  with  a one- 
meter  spectrograph.  The  top  spectrum  is  of  a krypton  lamp  used  for  wavelength 
calibration.  All  the  lines  in  the  spectrum  emitted  by  the  plasma  have  been  identi- 
fied as  originating  from  atomic  oxygen,  ionized  oxygen  or  hydrogen  atoms, 
indicating  a low  impurity  level. 
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Figure  l . Hybrid  schematic  of  the  50  k.I  o -pinch  used  as  a plasma  source. 


OXYGEN  SPECTRUM  FROM  50KJ,SL0W  0-PINCH 
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The  next  thing  done  was  to  study  the  time  evolution  of  an  atomic  oxygen  line 
and  an  ionized  oxygen  line.  This  was  done  using  a photomultiplier  tube  and  a one- 
meter  monochromator.  A plot  of  the  relative  intensities  of  these  two  species, 
versus  time  after  initiation  of  main  bank  discharge  is  shown  in  Figure  3.  It  is 
clear  that  although  the  emission  from  ionized  oxygen  dominates  the  spectrum  at 
early  times  the  atomic  oxygen  emission  will  be  the  only  emission  at  late  times. 

Since  these  preliminary  diagnostics  indicated  that  a slowly  decaying  neutral 
oxygen  plasma  had  been  generated  at  later  times,  more  detailed  diagnostics  were 
undertaken.  A Langmuir  (References  12  through  20)  probe  system  was  built  in 
order  to  measure  the  ion  density  and  electron  temperature  of  the  plasma  as  a 
function  of  space  and  time.  A schematic  of  the  Langmuir  probe  setup  is  shown 
in  Figure  4.  The  oscillator  is  a free-running  sawtooth  generator  with  adjustable 
frequency,  gain,  and  voltage  zero  level.  A free  running  oscillator  has  an  advan- 
tage over  a single  pulsed  system  in  that  the  probe  tends  to  be  cleaned  by  the 
large  currents  drawn  in  the  hot,  dense  early  discharge.  The  voltage  applied 
between  the  electrodes  is  measured  directly  while  the  current  flowing  is  found  by 
the  voltage  drop  across  a known  resistor.  Both  voltages  are  measured  differen- 
tially so  that  the  probe  is  free  to  float  at  the  plasma  potential. 

A picture  of  the  equal  area  double  probe  used  for  most  of  the  work  presented 
here  is  shown  in  Figure  5.  The  two  electrodes  of  the  probe  are  made  from 
tungsten  wire  0.040  inches  in  diameter  and  about  5 mm  long.  The  electrodes  are 
placed  ~2  mm  apart  which  is  the  approximate  spatial  resolution  of  the  probe. 

The  electrodes  are  mounted  on  a piece  of  bent  pyrex  tubing.  The  pyrex  tube  lies 
in  the  bottom  of  the  0 -pinch  tube  and  by  pulling  it  in  and  out  axial  variations  of  ion 
density  and  electron  temperature  can  be  measured.  Different  radial  positions 
can  be  investigated  by  rotating  the  pyrex  tube. 

A typical  double  probe  trace  is  shown  in  Figure  6.  The  upper  trace  is  the 
voltage  applied  by  the  ramp  generator  between  the  double  probe  electrodes.  The 
lower  trace  is  proportional  to  the  current  flowing  through  the  electrodes  of  the 
probe.  As  can  be  seen,  there  is  a saturation  of  the  current.  From  the  level  at 
which  the  current  saturates  1^  ^ one  can  determine  the  ion  density  N-  from 

i K eAp  \kTeJ 

where  e is  the  charge  on  an  electron,  A^  is  the  probe  area,  M-  is  the  ion  mass, 
k is  Boltzmann’s  constant,  Tp  is  the  electron  temperature,  and  K is  a constant 
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Figure  4.  A schematic  of  the  apparatus  used  to  make 
Langmuir  probe  measurements. 
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Figure  5.  The  equal  area  double  probe  shov  ing  the  two  tungsten  probe 
electrodes  and  a portion  of  the  long  pyrex  tube  upon  which 
it  is  mounted. 


chosen  to  be  0.55  as  in  Reference  16.  In  our  plasma  the  assumption  of  quasi- 
neutrality  should  be  good  and  thus  the  ion  density  should  equal  the  electron  density 
N^.  The  electron  temperature,  is  obtained  from  the  rate  of  change  of  current, 
I,  with  respect  to  voltage  V,  when  the  voltage  is  zero  as  given  by 


kTe 


(3) 


Three  radial  scans  and  many  axial  scans  were  done,  and  it  was  found  that  the 
temperature  was  not  a function  of  position  in  the  plasma  only  of  time.  The  elec- 
tron density  on  the  other  hand  was  found  to  vary  spatially  and  temporally.  Three 
radial  scans  at  the  coil  center.  44  cm  along  the  axis  from  the  coil  center,  and 
60.5  cm  along  the  axis  from  the  coil  center  are  shown  in  Figures  7 through  9. 
These  are  plots  of  the  ion  saturation  current  in  amps  (proportional  to  the  electron 
density)  versus  radial  position  from  the  tube  axis  at  various  times  after  initiation 
of  the  main  bank  discharge.  These  plots  show  a flat  electron  density  distribution 
over  the  central  4 cm  of  the  tube.  This  means  that  no  correction  need  be  made 
for  radial  inhomogeneities  as  long  as  the  light  collection  optics  is  stopped  down  to 
look  at  only  the  central  4 cm  of  the  plasma. 

A plot  of  ion  saturation  current  versus  axial  position  at  various  times  after 
initiation  of  main  bank  discharge  is  shown  in  Figure  10.  Zero  is  at  the  center  of 
the  single  turn  u -pinch  coil,  and  the  end  wall  is  at  the  end  of  the  quartz  tube.  In 
this  case,  there  is  an  inhomogeneity  and  a correction  mus*  be  made  for  it  since 
our  absolute  light  intensity  measurements  will  be  made  looking  axially  down  the 
tube.  Since  the  plasma  line  radiation  intensity  is  approximately  proportional  to 
the  square  of  the  electron  density  (Np) . the  easiest  thing  to  do  is  integrate  N* 
along  the  axis  of  the  tube.  One  can  relate  this  integral  to  the  square  of  the  elec- 
tron density  on  axis,  at  the  coil  center,  by  choosing  an  effective  length  ( ) 
which  when  multiplied  by  on  axis  at  the  coil  center  will  give  the  value  of  the 
integral.  A plot  of  versus  time  is  shown  in  Figure  11  . The  dashed  line  at 
~ 140  cm  is  the  true  length  of  the  tube.  Thus  all  the  homogeneity  corrections  are 
lumped  into  an  effective  length,  and  from  now  on  the  electron  density  on  axis  in 
the  center  of  the  coil  will  be  the  only  other  parameter  necessary  to  do  the  radia- 
tion transfer  for  our  optically  thin  plasma.  The  electron  density  and  temperature 
were  measured  on  axis  at  the  coil  center  as  a function  of  time.  Results  from 
these  measurements  will  be  compared  later  with  similar  results  from  Thomson 
scattering  measurements  in  Figures  14  and  15. 
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DENSITY  VS.  RADIAL  POSITION  AT  VARIOUS  TIMES 

AT  CENTER  OF  COIL 


REAR 


FRONT 


Figure  7.  Scan  of  ion  saturation  current  versus  radial  position  as  a 
function  of  time  after  initiation  of  main  bank  discharge. 
These  scans  were  taken  in  the  oxygen  plasma  in  the  mid- 
plane  of  the  single  turn  p -pinch  coil. 


DENSITY  VS.  RADIAL  POSH  ION  AT  VARIOUS  TIMES 
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Figure  8.  Scan  of  ion  saturation  current  versus  radial  position  as  a 
function  of  time  after  initiation  of  main  bank  discharge. 
These  radial  scans  were  taken  44  cm  down  the  tube  axis 
from  the  mid-plane  of  the  single  turn  ()  -pinch  coil. 
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60.5  cm  FROM  CENTER  OF  COIL 

I0"2, 

I L OAA  * . AAA  V W — ■ —A.  V 


sec/^x^- — * 


/ 250 


300 


350 

X--  X* 


/ 

/ 400, 

X* 

/ 45^ 
x/50p^ 

/s' 

/550 


x-«^  ,,x- 

' X — 


y , , x ’ 

'* -X-"~ 


- X x x. 


.X X- 


t 

\ 

\ 


-4  -3  -2  -I  0 

CM 


2 3 4 


Figure  9.  Scan  of  ion  saturation  current  versus  radial  position  as  a 
function  of  time  after  initiation  of  main  bank  discharge. 
These  radial  scans  were  taken  60.5  cm  down  the  tube 
axis  from  the  mid-plane  of  the  single  turn  P -pinch  coil. 
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The  apparatus  used  to  do  Thomson  scattering  (References  21  and  22)  is  shown 
in  Figure  12 . The  150  MW  ruby  laser  was  focused  on  axis  in  the  center  of  the 
coil.  The  scattered  light  was  viewed  at  a right  angle  to  the  incident  beam  through 
a slot  in  the  single-turn  main  coil.  The  scattered  light  was  detected  by  a photo- 
multiplier after  being  dispersed  by  a one-meter  monochromator.  Because  of  the 
reproducibility  of  the  0 -pinch,  we  were  able  to  scan  the  scattered  light  profile  on 

O 

a shot -to-shot  basis.  The  central  15  A of  the  scattered  light  profile  was  not 
scanned  since  the  stray  light  level  was  significant  in  this  region.  A typical 
scattered  light  profile  is  shown  in  Figure  13  This  is  a plot  of  relative  intensity 
versus  wavelength  taken  at  94  psec  after  initiation  of  main  bank  discharge.  The 
temperature  is  obtained  from  the  half  width  of  the  profile  while  the  electron 
density  can  be  obtained  from  the  absolute  total  intensity  of  the  scattered  light. 

The  absolute  total  intensity  of  the  scattered  light  was  found  by  calibrating  the 
photo- detect  ion  system  using  Rayleigh  scattering  from  molecular  oxygen. 

The  original  thought  behind  doing  both  Langmuir  probe  and  Thomson  scatter- 
ing diagnostics  was  to  use  the  Langmuir  probe  as  a relative  electron  density 
measure.  This  we  could  calibrate  at  early  times  in  the  discharge  by  doing 
Thomson  scattering  at  one  spatial  position  (on  axis  at  coil) , and  if  necessary  all 
the  Langmuir  probe  data  could  be  multiplied  by  some  constant  factor  to  make  it 
agree  with  Thomson  scattering.  As  can  be  seen  from  the  plot  of  electron  density 
versus  time  after  initiation  of  main  bank  discharge,  shown  in  Figure  14.  this 
procedure  was  unnecessary.  Similar  agreement  was  found  between  the  two 
diagnostic  techniques  when  the  electron  temperature  was  measured.  A plot  of 
electron  temperature  versus  time  after  initiation  of  main  bank  discharge  is  shown 
in  Figure  15. 

B.  Results 

A schematic  of  the  apparatus  used  to  make  absolute  intensity  measure- 
ments of  the  visible  bound-bound  radiation  of  oxygen  is  shown  in  Figure  16. 
Measurements  are  made  by  looking  axially  down  the  tube . The  aperture  stop  at 
the  first  lens  eliminates  the  need  for  radial  inhomogeneitv  corrections  as  dis- 
cussed in  the  diagnostics  section.  The  mechanical  shutter  prevents  light  emitted 
by  the  plasma  at  early  times  from  saturating  the  photomultiplier  tube  which  has 
a high  gain  for  detecting  weak  signals.  Calibration  of  the  detection  system  is 
done  in  situ  using  a standard  tungsten  lamp.  Absolute  intensity  measurements 
were  made  on  eight  visible  lines  of  the  neutral  oxygen  spectrum.  These  line 
wavelengths  along  with  the  transitions  (Reference  23)  from  which  they  are 
emitted  and  other  atomic  information  is  given  in  Table  1 . 
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Figure  13.  Atypical  scattered  light  profile  shown  fit  via  least 
squares  to  a Gaussian  ( solid  line) . The  error  bars 
represent  the  scatter  of  four  different  shots.  This 
particular  example  was  taken  94  pScC  after  initia- 
tion of  main  bank  discharge. 
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after  initiation  of  main  bank  discharge  taken  on  axis  in  the  mic 
plane  of  the  single  turn  n -pinch  coil.  The  agreement  between 
Thomson  scattering  and  Langmuir  probe  techniques  is  seen  to 
be  good. 
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Figure  15 . 


A plot  of  the  electron  temperature  of  the  oxygen 
plasma  versus  time  after  initiation  of  main  bank 
discharge  taken  on  axis  in  the  mid-plane  of  the 
single  turn  0 -pinch  coil. 
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Figure  16.  A schematic  of  the  apparatus  used  to  make  absolute 
line  intensity  measurements. 


TABLE  1 . OXYGEN  TRANSITIONS  STUDIED 


Multiplet 

Wavelength 

o 

(A) 


5330.0 

5436.1 

6157.3 


.4 


6455.0 

4368.3 
3947.29 

7773.4 


1 

Multiplet 

Transition 

Upper 

Level 

Probability 
1 06 ( sec-1 ) 

Degeneracy 

Ionization 
energy  (eV) 

"*D°  _ 

1 .97 

25 

0.554 

5 gO  5 p 

3.05 

5 

0.594 

^D°  — 

7.01 

25 

0.854 

3S°  - 3P 

6.20 

3 

0.914 

5 gO  5 p 

7.10 

5 

0.954 

3?  _ 3gO 

0.66 

9 

1 .254 

5p  _ 5S° 

0.33 

15 

1 .328 

5p  - 5S° 

34  .00 

15 

2 .874 

Figure  17  is  a plot  of  the  experimental  results  obtained  at  150  /xsec  after  the 

initiation  of  nain  bank  discharge.  The  number  density  divided  bv  the  degeneracy 

is  plotted  against  the  ionization  energies  of  the  various  levels.  At  these  condi- 

N « 2.8  x 10^  cm-^  and  T ^0.9  eV 
e e 

predicts  populations  of  the  energy  levels  should 


tions 


, the  code  developed  at  NRL 
be  in  Saha- Boltzmann  equilibrium 
with  the  free  electron  gas  (solid  line).  The  data,  however,  fits  very  well  the 
population  densities  expected  if  the  excited  levels  are  in  Boltzmann  equilibrium 
with  the  ground  state.  A possible  explanation  for  this  Discrepancy  is  that  since 
the  resonance  lines  are  optically  thick,  one  would  expect  the  excited  state  number 
densities  to  be  in  Boltzmann  equilibrium  with  the  ground  state  (References  24 
and  25) . The  reason  the  electron  number  density  is  not  in  equilibrium  with  the 
ground  state  is  that  there  is  a high  electron  loss  rate  out  the  ends  of  the  contain- 
ment vessel.  This  rate  is  over  a factor  of  ten  larger  than  the  electron  collisional- 
radiative  recombination  rate  and  tends  to  convert  electrons  plus  ions  into  ground 
state  atoms.  This  electron  sink  is  not  included  in  the  NRL  code.  The  error  in 
these  measurements  is  due  almost  entirely  to  the  uncertainty  of  the  absolute  value 
of  the  atomic  transition  probabilities  which  is  on  the  order  of  - 25  percent . This 
error  estimate  is  confirmed,  at  least  in  a relative  sense,  by  the  small  amount  of 
scatter  about  the  straight  dashed  line. 

The  problem  now  is  to  estimate  how  long  it  will  take  for  these  levels  to  relax 
from  a Boltzmann  equilibrium  with  the  ground  state  to  a collisional  radiative 
equilibrium  with  the  free  electron  gas  once  the  temperature  falls  sufficiently  to 
decouple  the  ground  state.  Knowing  the  optical  depth  and  the  excited  state  number 
densities,  a time  on  the  order  of  1 /xsec  was  estimated  to  be  sufficient  for  the 
triplet  system  whereas  the  relaxation  of  the  quintet  system  would  depend  on  the 
rate  of  spin  exchange  collisions.  As  can  be  seen  from  Figures  18  ;uid  19,  the 
relaxation  time  is  more  like  50  to  100  p sec  and  both  the  quintet  and  triplet  systems 
tend  to  relax  at  roughly  the  same  rate.  The  reason  for  the  slow  decay  rate  is  not 
understood,  but  the  equality  of  the  triplet  and  quintet  decay  rates  indicates  a large 
spin  exchange  rate.  It  also  is  apparent  from  the  figures  that  for  times  larger 
than  300  //.sec  after  initiation  of  main  bank  discharge  the  populations  of  the  excited 
states  will  be  determined  by  a collisional  radiative  model  and  will  not  be  radiative- 
ly  tied  to  the  ground  state. 

The  experimentally  determined  number  densities  of  eight  levels  as  a function 
of  ionization  energy  are  shown  for  300  //sec  after  initiation  of  main  bank  discharge 
in  Figure  20.  Figure  21  contains  a similar  plot  for  seven  levels  at  350  p sec  after 
initiation  of  main  bank  discharge.  The  uncertainties  of  the  individual  level  number 
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Figure  18.  Time  development  of  the  number  density 
of  the  atomic  oxygen  5^s  level. 
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Figure  20.  Comparison  of  experimental  number  densities 

of  oxygen  excited  states  divided  by  their  degenera- 
cies with  calculated  values  (0)  at  300  p sec  after 
initiation  of  main  bank  discharge. 
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density  measurements  due  to  measurement  error  alone  are  indicated  by  the  error 

bars  on  the  experimental  points.  The  error  bar  on  the  solid  line  in  each  figure 

represents  the  amount  of  cnange  which  can  be  obtained  in  the  collisional  model  by 

changing  the  electron  temperature  and  density  to  their  extreme  measured  values. 

The  number  densities  calculated  from  the  NRL  computer  program  are  also  shown 

in  the  figures.  The  program  uses  Seaton  type  cross-sections  for  collisional 

excitation  and  de-excitation  between  all  levels  with  dipole  allowed  transitions; 

i.e.,  n,  f i 1,  even  spin  exchange  collisions.  If  spin  exchange  collisions  are 

5 

ignored,  the  calculated  number  density  for  the  3 p level  is  almost  two  orders  of 

5 

magnitude  too  high.  This  indicates  that  the  3 s level  is  not  metastable  at  these 
conditions  and  can  be  collisional  depopulated  to  the  triplet  system  (References  26 
and  27) . It  also  should  be  noted  that  the  measured  number  densities  are  consider- 
ably lower  than  those  calculated  for  the  levels  which  have  an  ionization  energy 
between  0.9  and  1 .0  eV.  The  explanation  for  this  discrepancy  may  be  that  the 
code  does  not  take  into  account  the  influence  of  core  excited  states.  The  3s'  *D 
level  is  nearly  degenerate  with  the  levels  in  question  mid  can  be  depopulated 
rapidly  by  an  optically  thin  radiative  transition  to  the  2p  *D  level.  This  may  mean 
that  the  3s'  *D  level  acts  like  a sink  which  will  tend  to  depopulate  the  levels  with 
nearly  the  same  energy. 

IV.  NITROGEN  MEASUREMENTS 
A.  Diagnostics 

The  ambient  fill  gas  used  throughout  this  part  of  the  experiment  was  pure 
N2  at  a pressure  of  10  mTorr.  As  with  the  oxygen  work  discussed  in  the  previous 
section,  the  first  diagnostic  undertaken  in  nitrogen  was  to  tm\e  me  visible  spec- 
trum shown  in  Figure  22  . The  upper  spectrum  on  each  strip  is  neon  spectrum 
used  for  wavelength  calibration.  The  lower  spectrum  on  each  strip  shows  the 
neutral  and  ionized  lines  of  the  nitrogen  atom  clearly.  Besides  these,  however, 

-4- 

the  strong  bands  of  the  neutral  (Ng)  and  singly  ionized  molecule  ( Ng  ) can  also  be 
seen.  Of  course  nitrogen  is  well  known  for  having  strong  molecular  band  emis- 
sion in  the  visible  part  of  the  spectrum  (Reference  28).  The  only  lines  found 
which  were  not  associated  with  nitrogen  were  the  lines  of  the  hydrogen  Balmer 
series.  This  indicates  a high  degree  of  purity  in  the  plasma. 

•f-  4 

The  time  development  of  the  various  species  present  (N.  N . Ng . and  Ng  ) 
was  investigated.  It  was  found  that  the  N + , Ng.  and  N^  emission  occurred  only 
during  the  iirst  50  psec  of  the  discharge  and  after  that  the  emission  from  the 
neutral  atomic  nitrogen  lines  predominated.  Of  course  the  Ng  emission  bands 
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must  also  occur  very  late  in  the  afterglow,  but  this  emission  was  not  observed  on 
the  time  scale  of  this  experiment. 

As  in  the  previously  described  oxygen  work,  the  spatial  development  of  the 
nitrogen  plasma  as  a function  of  time  after  the  initiation  of  main  bank  discharge 
was  investigated  with  the  double  Langmuir  probe.  The  temperature  was  again 
found  to  be  constant  over  the  spatial  coordinates  and  vary  only  with  time.  The 
relative  ion  density  was  measured  as  a function  of  radial  position  for  various 
times  after  the  initiation  of  main  bank  discharge  at  the  mid-plane  of  the  single 
turn  main  coil.  A plot  of  the  results  of  this  measurement  is  shown  in  Figure  23. 
As  can  be  seen,  the  density  distribution  is  not  quite  as  homogeneous  as  in  the 
oxygen  case,  however,  it  is  still  felt  that  as  long  as  the  end-on  viewing  aperture 
is  held  down  to  4 cm  there  will  be  no  need  for  radial  homogeneity  corrections. 

The  axial  variation  of  ion  density  as  a function  of  time  after  main  bank  dis- 
charge was  also  investigated  with  the  free  running  double  Langmuir  probe. 

Figure  24  is  a plot  of  ion  saturation  current  (proportional  to  ion  density)  ns  a 
function  of  axial  position  for  various  times  after  main  bank  discharge.  As  can  be 
seen,  the  ion  density  is  very  nearly  constant  inside  the  main  coil  (0  to  40  cm) 
and  falls  off  somewhat  faster  at  the  later  times  outside  the  coil  than  the  oxygen 
plasma  did.  The  nitrogen  plasma,  as  did  the  oxygen  plasma,  also  displays  an 
ion  density  increase  at  the  end  of  the  main  coil  ' 250  p sec  after  the  initiation  of 
main  bank  discharge.  This  increase  is  probably  due  to  a reflected  shock  wave 
returning  from  the  end  wall  of  the  plasma  tube. 

An  effective  length,  (off.  was  also  calculated  for  this  plasma  in  the  same 
way  as  was  done  for  the  oxygen  plasma.  A plot  of  versus  time  after  the 
initiation  of  the  main  bank  discharge  is  shown  in  Figure  25.  As  can  be  seen,  the 
plasma  effective  length  tends  to  just  be  the  main  coil  length  of  80  cm  at  the  later 
times. 

The  electron  temperature  and  ion  density  were  measured  on  the  axis  of  the 
discharge  tube  in  the  mid-plane  of  the  main  turn  coil.  Plots  of  the  electron 
temperature  and  ion  density  as  functions  of  time  after  the  initiation  of  main  bank 
discharge  are  shown  in  Figures  26  and  27.  On  the  same  graphs  we  have  also 
plotted  a Thomson  scattering  measurement  of  the  electron  temperature  and 
density  done  100  p sec  after  the  initiation  of  main  bank  discharge.  The  Thomson 
scattering  value  of  the  electron  temperature  is  only  about  12  percent  lower  than 
the  Langmuir  probe  measured  value.  For  electron  density  Thomson  scattering 
gives  a value  about  20  percent  lower  than  found  from  the  Langmuir  probe  work. 
These  small  discrepancies  are  within  the  scatter  of  experimental  error  and  thus 
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Figure  23.  Scan  of  the  ion  saturation  current  versus  radial  position 
as  a function  of  time  after  the  initiation  of  main  bank 
discharge.  These  scans  were  taken  in  the  nitrogen 
plasma  in  the  mid-plane  of  the  single  turn  o -pinch  coil. 
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Figure  24.  A plot  of  ion  saturation  current  versus  axial  position 
for  the  nitrogen  plasma  at  various  times  after  the 
initiation  of  main  bank  discharge. 
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Figure  26.  A plot  of  the  electron  temperature  of  the  nitrogen 
plasma  versus  time  after  the  initiation  of  main 
bank  discharge  taken  on  axis  in  the  mid-plane  of 
the  single  turn  P-pinch  coil. 
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the  two  methods  again  show  very  good  agreement.  In  what  follows  the  electron 
temperature  and  density  will  be  used  as  measured  by  the  Langmuir  probe. 

B.  Results 

Absolute  intensity  measurements  of  eight  neutral  nitrogen  multiplets 
.vere  made  on  the  same  apparatus  used  <.o  take  the  previously  described  oxygen 
data.  Figure  16.  Table  2 gives  the  wavelengths,  transitions,  upper  state  degen- 
eracies and  ionization  energies,  and  the  atomic  transition  probabilities  of  the 
bound-bound  radiation  measured  (Reference  23) . The  number  density  divided  by 
the  degeneracy  of  the  upper  levels  of  these  eight  transitions  are  shown  at  150  p sec, 
200  nsec,  250  nsec,  300  / isec . and  350  psec  after  the  initiation  of  main  bank 
discharge  in  Figures  28  through  32. 

A o 4 

Note  should  first  be  taken  of  the  multiplet  produced  by  the  4p  - 3s  ( P - P) 
transition  shown  in  Figure  28  as  a low  point.  The  value  of  the  atomic  transition 
probability  for  this  transition  as  published  in  Wiese,  et  al.  , is  a factor  of  ten  too 
large  (Reference  29) . As  can  be  seen  in  the  figure,  a much  better  agreement  is 
obtained  when  the  number  density  of  the  upper  level  of  this  transition  is  multiplied 
by  ten.  In  Figures  29  through  31 , this  correction  has  been  carried  out,  but  not 
indicated.  The  uncertainty  in  the  transition  probabilities  for  these  eight  nitrogen 
transitions  is  about  50  percent,  as  opposed  to  ^25  percent  for  the  oxygen  transi- 
tions measured,  and  this  larger  uncertainty  is  indicated  by  the  larger  scatter 
about  the  theoretical  curve  in  Figure  28. 

The  next  thing  to  notice  about  the  150  nsec  measurements  is  that,  as  in  the 
case  of  the  oxygen  studies,  the  ground  state  is  over-populated.  However,  in  this 
case,  the  ground  state  is  over-populated  by  about  a factor  of  500  instead  of  the 
factor  of  3 found  for  oxygen.  This  is  due  to  the  higher  electron  temperature  ;uid 
density  of  the  nitrogen  plasma.  When  one  uses  the  ground  state  over-population 
factor  along  with  the  measured  temperature  to  interpolate  between  the  curves 
published  by  Park  (Reference  25)  from  his  radiative-collisional  non-equilibrium 
calculation  on  nitrogen  one  obtains  the  dashed  line  on  Figure  28.  The  good 
agreement  is  rather  fortuitous  since  a small  error  in  the  plasma  temperature 
will  make  a large  difference  in  the  calculated  over-population  of  the  ground 
state.  When  a reasonable  error  is  allowed  in  the  temperature  measurements 
made  at  200  p sec  after  the  initiation  of  main  bank  discharge,  one  finds  that  the 
data  at  this  time  is  also  consistent  with  the  calculations  done  by  Park. 

A problem  arises  when  the  measurements  250  psec.  300  p sec.  and  350  p sec 
after  main  bank  discharge  are  considered.  The  ground  state  is  several  orders  of 
magnitude  under-populated  at  those  times  relative  to  the  equilibrium  population 
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TABLE  2.  NITROGEN  TRANSITIONS  STUDIED 


Multiplet 

Transition 

Upper  Level 

Wavelength 

0 

(A) 

Transition 

Multiplet 

Probability 
10®( sec"^ ) 

Degeneracy 

Ionization 
energy  (eV) 

5201 .80  | 

5197.80  J 

2p25d  - 2p23p 

2P  - 2S° 

2.30 

6 

0.547 

6646.51  ) 
6644.96  j 

2p25s  — 2p23p 
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Figure  28.  A plot  of  tlio  experimental  number  densities  of 
nitrogen  excited  states  divided  by  their  degen- 
eracies versus  their  ionization  energies  at  150 
psee  after  the  initiation  of  main  bank  discharge 
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Figure  29.  A plot  of  the  experimental  number  densities  of  nitrogen 
excited  states  divided  by  their  degeneracies  versus 
their  ionization  energies  at  200  psec  after  the  initiation 
of  main  bank  discharge. 
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Figure  30.  A plot  of  the  experimental  number  densities  of  nitrogen 
excited  states  divided  by  their  degeneracies  versus 
their  ionization  energies  at  250  jisec  after  the  initiation 
of  main  bank  discharge. 
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Figure  31.  A plot  of  the  experimental  number  densities  of  nitrogen 
excited  states  divided  by  their  degeneracies  versus 
their  ionization  energies  at  300/xsec  after  the  initiation 
of  main  bank  discharge. 
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Figure  32.  A plot  of  the  experimental  number  densities  of  nitrogen 
excited  states  divided  by  their  degeneracies  versus 
their  ionization  energies  at  350  p sec  after  the  initiation 
of  main  bank  discharge . 


and  the  levels  measured  should  fall  below  the  Saha-Boltzmann  line  drawn  on  the 
plots.  However,  as  can  easily  be  seen,  the  levels  are  over-populated  and  this 
trend  gets  larger  at  later  times.  This  effect  seems  real  since  the  discrepancy 
lies  outside  the  indicated  errors  on  the  measurements,  and  even  invoking  the  50 
percent  uncertainty  on  the  individual  transition  probabilities  does  not  bring  agree- 
ment. This  observation  indicates  that  some  reaction  is  occurring  which  supplies 
excited  state  nitrogen  atoms  at  a higher  rate  than  collisional-radiative  recombina- 
tion. One  possible  energetic  reaction  is  N + N + N - N2  + N*  where  N*  indicates 
an  excited  nitrogen  atom.  In  order  for  this  reaction  to  compete  with  electron 

radiative-collisional  recombination,  it  must  have  a rate  coefficient  on  the  order 
- 1 7 3 

of  3 x 10  cm  sec.  This  is  a bit  larger  than  presently  accepted  but  other 
energetically  possible  reactions  fall  much  farther  short  of  the  needed  rate. 

Another  possibility  that  was  considered  was  that  the  Langmuir  probe  intro- 
duced impurity  oxygen  into  the  nitrogen  plasma  thereby  lowering  the  electron 
density  at  the  late  times  through  the  reactions 

N + + 02  - N + C>2 

and 


°2  + 0 ~ °2  * 

Then  when  measurements  of  the  absolute  line  intensity  were  made,  without  the 
probe  being  present,  the  electron  density  was  higher.  This  explanation  can  only 
be  true  if  radiative-collisional  recombination  through  the  impurity  channel  con- 
trols the  electron  loss  rate  in  the  plasma.  The  electron  loss  rate  for  this  plasma 
is  completely  dominated  (by  several  orders  of  magnitude)  by  the  flow  of  electrons 
and  ions  out  the  open  ends  of  the  P-pinch  tube.  This  means  that  even  the  addition 
of  a large  amount  of  impurity  would  not  affect  the  electron  density  of  the  plasma. 
V.  SUMMARY 

14  -3 

The  early  time,  (Np  - 10  cm  Tp  ^ 1 eV)  oxygen  and  nitrogen  results 
are  well  understood  when  tiie  over-population  of  the  ground  state  is  taken  into 
consideration.  The  early  time  nitrogen  results  agree  very  well  with  the  calcula- 
tions done  by  Park,  but  the  NRL  radiative-collisional  recombination  code  does 
not  seem  to  giv^  the  correct  results  in  this  limit  for  oxygen. 

At  late  times,  (Np  ~ 10^  cm  ^ . Tp  - 0,4  eV)  the  NRL  code  gives  good 
agreement  with  the  oxygen  experiment  if  large  collisioaal  spin  exchange  cross- 
sections  are  used.  However,  the  NRL  code  does  not  reproduce  the  experimentally 
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observed  under-population  of  the  levels  with  ionization  eneigies  of  ~0.9  eV  since 
it  does  not  take  into  account  the  depopulizing  effect  of  the  3s’  level.  The  late 
time  behavior  of  the  nitrogen  plasma  is  not  well  understood  at  this  time  (Refer- 
ence 30) . However,  since  the  discrepancy  is  larger  than  can  be  accounted  for  by 
probable  er.perimental  inaccuracies,  we  believe  it  to  be  real. 
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